Field dependent equilibrium energetic distribution of localized charge carriers in 
disordered semiconductors at low temperatures. 
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The hopping transport of charge carries in 3-dimensional 
disordered semiconductors is analyzed in this work. It is as- 
sumed that the density of localized states is a decreasing func- 
tion of energy and contribution of thermoactivated jumps into 
the transport is negligible. The possibility of equilibrium en- 
ergetic distribution for this case is shown. This distribution 
is characterized by field-dependent effective temperature. 

Processes of carrier transport in amorphous semicon- 
ductors with conduction via dclocalized states have been 
successfully described in terms of the multiple-trapping 
model [1-3]. This model implies that two carrier frac- 
tions exist in the material, occupying at any given in- 
stant of time delocalized and localized states. How- 
ever, in the materials with all carriers in localized states, 
e.g., in most polymers, charge transport is due to car- 
rier hopping. Many interesting and important features of 
the hopping transport in disordered semiconductors have 
been observed experimentally in fairly high electric fields 
[4-6]. It has been shown both analytically [7-10] and by 
Monte Carlo simulations [11, 12] that the energy distri- 
bution of localized carriers and the carrier mobility in a 
disordered semiconductor can be described both by the 
ordinary temperature T and by some effective tempera- 
ture Tp(F) dependent on the magnitude of the applied 
external electric field F. Analytical calculations [7, 9] of 
this effective temperature yield 
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while numerical simulation yields Tp = 0.67 • eF/jk (see 
[11]) or Tf = (0.69±0.03)-ei^/7/c (see [12]), where e is the 
elementary charge, I/7 the localization radius, and k the 
Boltzmann constant. The equilibrium energy distribu- 
tion of injected carriers in an amorphous semiconductor 
subjected to an external electric field at low tempera- 
tures, when the contribution of thermally activated hops 
to the transport process is negligible, has been found [9] 
for the ID case, which, in fact, corresponds to high elec- 
tric fields. In this paper, we demonstrate that the equi- 
librium energy distribution of localized carriers is also 
possible in the 3D case at lo w temperatures in an exter- 
nal electric field. We show that this distribution can be 
described by the Boltzmann exponential with an effective 
temperature given by (1). 

As an equation describing the kinetics of the hopping 



transport, we use the well-known balance equation 
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written for the case of low occupancy of localized states, 
fi ^ I. Here, t is the time, fi the mean occupation 
number of the i-th state and Vij the probability of tran- 
sition from the state i to the state j . Let us assume 
that in the quasiequilibrium transport regime, when the 
density of states is time-independent, a great number 
of carriers can a void unlikely jumps that take a very 
long time tfi ^ h/AE, i.e., jumps that do not con- 
tribute significantly to the transport process [13]. Here, 
h is Planck's constant and AE the characteristic energy 
change in carrier transition from one localized state to an- 
other. This assumption makes it possible to take advan- 
tage of the concept of the distribution function / (r, E, t) 
averaged over continuous variables: energy E ("deeper" 
states have higher energies E) and radius vector r. Let 
us use the Miller- Abrahams expression for the transition 
rate 
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where £ = E — eFAr — E' , H{£) is the Heaviside func- 
tion, vq the attempt-to-jump frequency, and Ar = r — r'. 
Then, we use equation derive the following kinetic equa- 
tion for the distribution function /, 
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by the condition 



dv f{r,E,t)dE=l, 



where g{E) is the energy distribution of locaHzed states 
(DOS) normalized by the condition 



g{E)dE = 1. 



(6) 



The zero energy E = corresponds to the maximum of 
the function g{E). Equation (4) is written on the as- 
sumption that the position and energy of locahzed states 
are uncorrelated, which corresponds to the case of com- 
pletely disordered materials. 

At low temperatures, when the contribution of ther- 
mally activated hops can be neglected subject to the con- 
dition of hypothetical equilibrium, kinetic equation (4) 
for the equilibrium distribution function averaged over 
the space coordinates 
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takes the form 
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We use equation (8) and formulas (5) and (6) to derive 
the following integral equation: 
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Here, we introduced the integrated operator 
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and the dimensionless variables e = E/kTp and e' = 
E'/kTp, where Tp is the effective temperature defined by 
(1). The appearance of this effective temperature is as- 
sociated with the fact that, during transport, the electric 
field transfers carriers into states with a higher energy, 
similarly to the thermodynamic effects described by the 
ordinary temperature T. 



This equation can be solved by the iteration method; 
as the first approximation, we use the density of localized 
states 



aE)=9{E). 
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(5) Substituting (11) into (9), we obtain 
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Here, the quantity EmiTp) depends on the particular 
form of the DOS function. In particular, in the case of a 
Gaussian distribution 

g{E)=goex.p{-[E/Eo]''), 

we have E^ (Tp) = Eq/2Tp, while in the case of an 
exponential distribution 



g{E)=goexp{-\E/Eo\) 

we obtain Em (Tp) = E^/{Tp - Eq). 
Hence, 
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Thus, to a first approximation, the equilibrium energy 

distribution is described by the Boltzmann fmiction. To 
calculate the distribution function feq{E) more precisely, 
we substitute, in accordance with (13), the function 
g{E)exp{E / kTp) into the right hand side of equation 
(9). The resulting solution feq (E) has the same asymp- 
totic form as the function feq^E), since the asymptotic 
form (12) of the function A {e, e') g{e') appearing in equa- 
tion (9) is not changed. Thus, if the DOS fimction g{E) 
is sufficiently "shallow", so that the equilibrium can be 
established, 
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then the asymptotic form of the occupancy of lo- 
calized states (f{E) = f{E)/g{E) is the exponential 
exp{E/kTp). The quantity E-miTp) appearing in equa- 
tions (12) and (13) is the characteristic energy of the peak 
in the distribution function feq{E). 

Thus, we showed the following. First, in the 3D case, 
as also in the ID one [9], the spatially uniform equilib- 
rium energy distribution of localized charge carriers may 
occur even in the limiting case T — > 0, when the con- 
tribution from thermally activated hops to the kinetic 
equation (4) can be completely neglected. Second, this 
distribution can be approximated by the Boltzmann dis- 
tribution function with the effective temperature (1) in- 
stead of the ordinary temperature. It should be noted 
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that this effective temperature coincides with previous 
analytical estimations [9, 10]. 
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